Anomalies and non-log-normal tails in one-dimensional localization with power-law disorder 
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Within a general framework, we discuss the wave function statistics in the Lloyd model of Anderson localiza- 
tion on a one-dimensional lattice with a Cauchy distribution for the random on-site potential. We demonstrate 
that already in leading order in the disorder strength, there exists a hierarchy of anomalies in the probability 
distributions of the wave function, the conductance, and the local density of states, for every energy which cor- 
responds to a rational ratio of wave length to the lattice constant. We also show that these distribution functions 
do have power-law rather then log-normal tails and do not display universal single-parameter scaling. These 
peculiarities persist in any model with power-law tails of the disorder distribution function. 
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Wave function localization in one spatial dimension has 
attracted enormous attention since the pioneering work by 
Anderson. 1 Much of our present understanding of this phe- 
nomenon is based on the original Anderson model taken on a 
one-dimensional lattice 



-t(#„+i + ¥„_i) + V n q„ = E$ n , 



(1) 



with a white-noise disorder (V n V m ) cx S nm , (V n ) = 0, and 
fixed hopping element t. The potential V n at each site takes 
real values according to a probability density Py(V). Most 
of the theoretical investigations of Anderson localization as- 
sume a finite variance var V = 2D < oo, and then con- 
sider the weak-disorder limit D <C t 2 . This condition is a 
prerequisite for single-parameter scaling, 2 in which the prod- 
uct L£~ x of system length L and inverse localization length 
= — ]iin n _ >0O — ln|^ n | is the only free parameter in 
the universal distribution function of the Lyapunov exponents 
a = — — In \ ^ n \ f° r finite n. This carries over to universal 
distribution functions of the dimensionless conductance g and 
the local density of states v. For L/£ 3> 1 these distribu- 
tion functions follow log-normal laws — the paradigm for 
large fluctuations. On the other hand, in many realistic appli- 
cations the distribution function Py(V) displays power-law 
tails, with D = oo, such that in a sense disorder never really 
is weak. The most prominent example is the localization of 
wave functions in the momentum space of the kicked rotator. 3 
This dynamical problem has been mapped onto the Anderson 
model in the seminal works I3l4l with an effectively random 
Cauchy-distributed potential 



Pv, Cauchy (V) = — Im — - 

7T V - iS 



8>0. 



(2) 



Equation Q with the disorder distribution function given by 
Eq. Q has been proposed for the first time by Lloyd. 5 In this 
model the localization length can be computed analytically for 
arbitrary S££ and the variance of the Lyapunov exponents a 
has been analyzed very recently in Ref.|8| 

A beautiful experimental realization of the kicked rotator 
is the dynamics of atoms driven by a regular train of laser 
pulses. 910 In these experiments, the probability distribution 
function in momentum space is seen to relax from an initial 
Gaussian into an exponential profile, demonstrating the ab- 
sence of diffusion in momentum direction. However, since 



the wave function statistics in localization is a prototypical 
example of large fluctuations, the localization length of the 
eigenstates can be quantitatively inferred from moments of 
the wave function (like the measured mean probability) only 
if the full shape of the distribution function is known, includ- 
ing its tails. The impact of the tails of the probability density 
Pv (V) on the statistics of the wave functions has been already 
mentioned by Halperin, 1 1 but has not been analyzed, let alone 
sufficiently appreciated, in the literature. 

In this paper we address the fluctuations of the wave func- 
tion and related quantities for distributions Py(V) with 
power-law tails, by going beyond the mean of a and its 
variance, studied so far^fi^ We first set out a general frame- 
work for arbitrary Pv(V), which then is applied to the Lloyd 
model with Py (V) given by Eq. The fluctuations turn out 
to be highly non-universal, with an anomalous energy depen- 
dence [reflecting also the spatial discreetness of the Anderson 
model Q] and non-log-normal tails that strongly affect the 
behavior of the moments of g, and v even in the weak- 
disorder limit S <C t, and even for low fractional orders of the 
moments. These characteristics of the wave function statis- 
tics are in striking contrast to the universality for models with 
D <C i 2 and single-parameter scaling. 

The central quantity of interest in our calculation is the gen- 
erating function /i(A) of the cumulants of In ^. As pointed out 
by Borland 12 and Thouless, 13 instead of solving Eq. Q as a 
boundary-value problem it suffices to investigate the specific 
solution <!>„ of the initial-value problem <I>o = a, $i = b. At 
large distances n 3> 1 this solution exponentially increases 
as $„ ~ exp(cm) for almost all values of a and b, which 
statistically is equivalent to the inverse wave-function decay 
^„ ~ Q^ 1 in the original problem. The cumulant-generating 
function 



/i(A) = lim — In 

n — *oo Ji 



fc=i 



c fe X k 



(3) 



accounts for the details of convergence of the Lyapunov expo- 
nent a to its mean value c\ = The coefficients Ck with 
k > 2 are numerical constants which characterize the devi- 
ations of a from c\. As follows from Eq. the cumulants 
of a vanish according to the law dictated by the generalized 



2 



central-limit theorem 



Ckn 1 h for 7i > £. 



(4) 



The coefficients Ck do not depend on the initial conditions for 
4>„ and capture the universal information about the fluctua- 
tions of many essential quantities in the localized regime, such 
as the dimensionless conductance g of a finite sample of the 
length L and the local density of states 1 ^ v in a semi-infinite 
wire at a distance L to an open boundary. Both quantities do 
strongly fluctuate even in the localized regime L>(, unlike 
the Lyapunov exponent. However, the fluctuations of their 
logarithms can be expressed through the same coefficients Ck 
by 

lim L(((-\ng) k ))= lim ±«(- W) fe )) = 2 k c k . (5) 

L — >oo Li L — >oo J_i 

Each coefficient may be used to define a length scale = 
c^ 1 , and the question for single-parameter scaling can be 
posed as whether these length scales are independent quan- 
tities or notiA£ 

We now set out a general approach to calculate the gener- 
ating function /.i(A) and the coefficients for arbitrary form 
of Pv(V). We built up on the formalism previously used to 
calculate the inverse localization length = c\ The 
Anderson model Q can be written in terms of new variables 
z n = * n+ i/* n , r n = In in the simple way 



%n 1/ %n- 



•ln|z n _i|, (6) 



where v n = (V n — E) jt. We seek the specific solution of the 
initial-value problem zq = b/a, ro = In \a\. Iterating the map 
l|6} we observe that z n and r n take real values z and r with a 
probability density P n (z, r), which obeys 



P n (z,r) = / F(v)dv // P n ^(z',r')dz'dr' 

J — oo JJ—oo 

x S(z-v + l/z')S(r-r' -ln\z'\), 

with F(v) = tP\/{vt + E) the probability density of v. 
It is convenient to introduce the function 



(7) 



h n (z,X) = \z\ x / e rX P n ( Zl r)dr 

J —oo 

and rewrite Eq. Q as 

h n (z, A) = |z| A / F(z + 1/z') K-itf, A) dz'. 



According to Eqs. and (jSJl we have 



M (A) 



lim — In 

n — >oo fi 



dz h n —i(z, A) 



(8) 



(9) 



(10) 



If the function /i(A) exists, the solution to Eq. at large n 
must fulfill the relation h n (z, A) = h n ~i(z, A). Thus 
Eq. (|9jl is transformed into the functional eigenvalue problem 

e nW-^\z\h( z ,\) =F[h](z,\), (11a) 

p OO 

F[h](z,\)= / F(z + l/z')h{z',X)dz'. (lib) 



This is the central general equation of this paper. In any prac- 
tical case it has to be solved perturbatively in A. We expand 
the function h(z, A) in a series 



A 



h{z,X) = -^i h k( z ) 



k=0 



(12) 



and introduce the notation c±(z) = Ci — In \z\. Equation dl It 
is transformed into the following set of equations 

F[ho] -ho=0, (13a) 

!F[h{\ — h\ = ci/io, (13b) 

F[h 2 \-h2 = (cl + c 2 )ho+2c 1 h 1 , (13c) 
F[hJ[ -h 3 = (cf + 3cic 2 + c 3 )h 

+3(5? + c 2 )hi + 3ci/i 2 , etc. (13d) 

Equation dl 3ai delivers the stationary distribution function of 
z and has been used beforai^iii^ to calculate the localization 
length £ = c^ 1 from 



ci 



/io(z) In |z c?2;. 



(14) 



So far we have shown that Eq. < !13al can be considered as just 
the first member of a hierarchy of equations that determine the 
complete wave function statistics for finite n 3> £. 

The integrals <fe of the left-hand sides of Eqs. (1 1 31 
equal zero. Equation dl4> indeed can be derived by integrat- 
ing both sides of Eq. (II 3bi along the real axis. Once the dis- 
tribution function ho(z) and the mean Lyapunov exponent c\ 
are known (as is analytically the case in the Lloyd model), one 
can construct the solution to Eq. dl3b> iteratively by 19 

hx(z) =K[c l ho](z), (15a) 

/OO 
K(z,z')f(z')dz'. (15b) 
-CO 

The second coefficient C2 in the cumulant expansion is readily 
found by integrating Eq. < ll 3cl along the real axis, 

/oo 
[h (z) (In \z\ - ci) + 2h 1 (z)](\n \z\ - Cl )dz. (16) 
-oo 

This procedure can be repeated recursively to calculate all 
other coefficients c^. 

For the Lloyd model, the scheme developed above allows 
us to obtain the coefficients Cfc analytically. The stationary 
distribution function of the variable z is readily found by the 
exact solution of the integral equation dl 3a> .^i^ 

h (z) = -Im— !— , s + - = E+tS ) lms>0. (17) 

IT Z — S S t 

The integral dl4> yields the well-known result 



ci = In Isl 



(18) 
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FIG. 1 : The ratio C3 / ci is plotted according to the analytical result 
12 1 1 for the Lloyd model at the disorder strength S = O.Olt. The ratio 
is never small inside the band and reveals anomalies at energies E = 
—2tcos(np/q) with p and q integers. The corresponding rational 
number p/qis indicated in the figure. The size of the anomaly only 
depends on the value of the denominator q. 



The kernel function K (z, z') can be obtained by iterative ap- 
plication of the operator J-[h], 



K{z,z') 



IlmVf 1— 



1 



Pn{z) 



(s n - s- n ) - z(s 



n+1 



z-r„ 

,-(«+!)' 



p n (±oo), 



,-(n-l) 



z(s T 



(19a) 

(19b) 
(19c) 



where the second term in the parenthesis on the right hand 
side of Eq. ( I19ai is added to provide a better convergence of 
the intermediate expressions^ 

Applying the result dl9l to Eqs. dl5> and d!6t and perform- 
ing the summation one recovers the result of Ref . [H 

Re [Li 2 (s- 2 ) - Li 2 (|s|- 2 )] + arg(s)(^ - arg(s)) 



(>2 



In |s| 2 [ln(|s| 2 - 1) -ln|s 2 - 1| 



(20) 



where Li„(z) = X^i z k /k n is the polylogarithmic func- 
tion. As has been shown in Ref. |H in the limit <5 — > the 
ratio c 2 jc\ equals 2 (not 1 as for conventional weak disorder) 
inside the band (it vanishes outside the band). This energy- 
insensitivity has encouraged the authors of Ref. to conclude 
that single-parameter scaling is fulfilled, and to introduce a 
novel criterion for single-parameter scaling. As we will dis- 
cuss now, these findings do not carry over to the fluctuations 
beyond the variance, characterized by c,t with k > 3. 
The coefficient C3 can be found from Eq. dl 3dl as 



c 3 = 3 



dz In 

00 



£(*) 



z ■ 

I S I 

dz In — 

!~W-i ' Z I 



— — im 

n— 1 



1/s 
In 



c 2 - In 2 



-2 In 



1 



-Pn(s*) 



1/8 
hi 



h (z) 

£(z), (21a) 
(21b) 





D ilA t Yli < e < t 


D/t < \E\ -C t 


e, \E\ ~ t 


Cl 


C2 


D/(Ue) 


D/{tt 2 ) 


D/(it 2 - E 2 ) 


C3 


C4 


33D 3 /(128eH 2 ) 


9£>7(32 E 2 t 4 ) 


oc D 3 /t 6 


Co 


C6 


5175D 5 /(2048£ 7 t 3 ) 


135D 5 /(128 E 4 t 6 ) 


<xD 5 /t 10 



TABLE I: The leading asymptotic values of the coefficients in 
the case of weak Gaussian disorder D <C t 2 upon the deviation 
e = 2t — \E\ from the band edge (second column) or from the band 
center (third column). These results are obtained from the saddle- 
point analysis of Eq. Jl H . 22,23 The last column represents the generic 



values inside the band. 
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FIG. 2: Distribution function P\ n g (In g) obtained from the Anderson 
model Q with Cauchy disorder (full circles), box disorder (open cir- 
cles), and Gaussian disorder (open squares). Parameters are in the 
localized regime L/£ = 6, with localization length £ = 200 lattice 
constants, and E = Q.2t. In this semilogarithmic plot, the lognormal 
distribution function for conventional weak (box or Gaussian) disor- 
der maps to an inverted parabola, while the straight-line asymptotics 
found for Cauchy disorder correspond to a power-law tail in P a {g). 



where the function p^ n 1 (z) stands for the inverse of p n (z). 
The ratio c^/ci is plotted at Fig. Q versus the energy. The 
plot clearly displays a sequence of sharp dips, which appear 
exactly at energies E = —2t coslnp / q) where p and q are 
integer, and become more narrow in the limit 6 — > 0. The 
anomaly in the band center is the biggest one and reaches 
about 3% of the absolute value of the ratio c^/ci in the limit 
5 — > 0. The existence of such anomalies for the inverse local- 
ization length = c\ has been pointed out by Lambert, 20 
but for this quantity they only show up in higher orders of the 
expansion in the disorder strength, with exception of the band 
edge \E\ = 2 and the band center E = 0pii22i2iFor conven- 
tional weak disorder with D <C t 2 the other anomalies should 
be seen in the higher coefficients with k > 3. However, 
those cumulants are themselves suppressed by orders of D jt 2 
(see Table [j}, again with the exception of the band edge and 
the band center, where they are of the same order as C\ and 

„^ 15.22.23.24 

In striking contrast, in the Lloyd model the coefficients 
Cfe increase very rapidly with increasing index k. In the 
limit S — > we indeed observe c 2 /ci = 2, c^/c\ = 5, 
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C4,/c% « 20, cg/ci w 100 for p/q irrational. The analysis 
of Eqs. dl 31 for the Lloyd model demonstrates that the gen- 
erating function /i(A) exists only for A < A c , where the con- 
vergence radius A c < 1, which implies a factorial growth of 
the ratios c k jc\ for large k. Such a behavior is consistent 
with a power-law tail in the conductance distribution function 
Pg(d) ~ g~( 2 ~ A c)/ 2 for g — > 0. For a general power law 
Pv{V) oc for \V\ — > oo, A c < (3 - 1 must be expected 

to depend on /3, implying that the precise form of the tail in 
P g {g) is not universal. The power-law tail in P g (g) is certi- 
fied by the numerical result for the conductance distribution 
function for Cauchy disorder, shown for L/£ = 6 in Fig.|2] 
The probability to find a vanishing conductance in the Lloyd 
model is strongly enhanced as compared to the case of con- 
ventional weak disorder, which displays log-normal tails. In 
this case, the generating function /x(A) is well-defined for all 
A in the whole energy range. Moreover, far from the band 
edges (e = It - \E\ > D 2/3 t~ 1/3 ) and far from the band 
center (E ^> D/t), it acquires a universal parabolic form 
fi(X) = £ _1 (A+A 2 /2), since c x = c 2 = D/(4t 2 -E 2 ), while 



all other coefficients can be disregarded in the limit D <C t 2 
(see Table[I}. 

In conclusion, we have studied analytically the statistics 
of localized wave functions in the Lloyd model, which is 
frequently used to analyze dynamical localization. We have 
found that the distribution functions of the conductance g and 
of the local density of states v do not have a log-normal form. 
Moreover, even in the limit of vanishing disorder these dis- 
tribution functions reveal sharp anomalies at energies E — 
—2tcos(irp/q), with p/q a rational number. These specific 
features can be attributed to the power-law decay of the disor- 
der distribution function. They sensitively affect the moments 
(including fractional moments) of g and v and demonstrate 
that for such distribution functions the wave-function statis- 
tics are highly non-universal. It would be striking to see sim- 
ilar effects for the prelocalized states in the diffusive regime 
L < £ of multichannel systems^ 

We thank B. L. Altshuler and S. Fishman for illuminating 
discussions. 
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